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Abstract Double dot exhange only qubit represents a promising compromise 
between high speed and simple fabrication in solid state implementations. 
A couple of interacting double dot exhange only qubits, each composed by 
three electrons distributed in a double quantum dot, is exploited to realize 
Controlled-NOT (CNOT) operations. The effective Hamiltonian model of the 
composite system is expressed by only exchange interactions between pairs 
of spins. Consequently the evolution operator has a simple form and repre¬ 
sents the starting point for the research of sequences of operations that realize 
CNOT gates. Two different geometrical configurations of the pair are con¬ 
sidered and a numerical mixed simplex and genetic algorithm is used. We 
compare the non physical case in which all the interactions are controllable 
from the external and the realistic condition in which intra-dot interactions 
are fixed by the geometry of the system. In the latter case, we find the CNOT 
sequences for both the geometrical configurations and we considered a qubit 
system where electrons are electrostatically confined in two quantum dots in 
a silicon nanowire. The effects of the geometrical sizes of the nanowire and of 
the gates on the fundamental parameters controlling the qubit are studied by 
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exploiting a Spin Density Functional Theory based simulator. Consequently, 
CNOT gate performances are evaluated. 

PACS 03.67.Lx • 73.21.La • 75.10.Jm 


1 Introduction 

Several solid state physical implementations of quantum computing have been 
proposed both from experimental P5M1I51IM7] and theoretical [RIIMTU] points 
of view. In solid state physics many approaches have been used ranging from 
quantum dots, to donor-atom nuclear spins or electron spins. A characteris¬ 
tic common to all the scenarios is a tunable Heisenberg exchange interaction 
between spins to generate two-qubit quantum gate operations. 

In the last decade electron spins confined in quantum dots have become an 
attractive basis for quantum computing because of their relatively long coher¬ 
ence times and potential for scaling mums- In the simplest proposal, 
single spins form the logical basis, with single qubit operations via spin reso¬ 
nance mm- An alternative scheme, with logical basis realized with singlet 
and triplet states of two spins requires inhomogeneous static magnetic field for 
full single-qubit control Huang- Using three spins to represent each qubit 
removes the need for an inhomogeneous field m- Exchange interactions be¬ 
tween adjacent spins suffice for all one- and two-qubit operations m- Here 
we study a pair of double dot exchange only qubits in which three electrons 
are distributed in a double quantum dot WMl- The logical states are de¬ 
fined adopting linear combination of singlet-triplet states of a pair of electrons 
in one dot with single electron states of the electron in the other dot. The 
Schrieffer-Wolff method has already been applied to a single qubit {221. We 
employ the method to treat the two qubits case in order to express the Hamil¬ 
tonian in terms of effective Heisenberg exchange and to exploit it to generate 
CNOT gates. The effective Hamiltonian explicitly express the dynamical be¬ 
havior of the system and determines the conditions to realize gate operations. 
The most efficient sequences of interactions realizing CNOT gates have been 
consequently calculated in two different geometrical configurations. A numer¬ 
ical calculation based on a simplex method mixed with a genetic algorithm 
has been used. Device designs featuring different geometrical sizes are com¬ 
pared by exploiting a Spin Density Function Theory (SDFT)-based simulator, 
presenting design rules for fast CNOT gates. 

The paper is organized as follows. Section 2 is devoted to the derivation 
of the effective Hamiltonians for the pair of interacting qubits in the two con¬ 
figurations under study. Section 3 contains the main results on the sequences 
that realize CNOT gates in the case in which all the interactions are control¬ 
lable from the external and in the realistic situation in which the intra-dot 
interactions are fixed by the geometry of the system. Finally in Section 4 an 
analysis on the gate performance in a silicon nanowire double quantum dot is 
presented. 
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2 Effective Hamiltonian for two interacting qubits 

In this Section we derive the Schrieffer-Wolff effective Hamiltonian to de¬ 
scribe in a simple and intuitive way a couple of interacting qubits by com¬ 
bining a Hubbard-like model with a projector operator method. As a result, 
the Hubbard-like Hamiltonian is transformed into an equivalent expression in 
terms of the exchange coupling interactions between pairs of electrons. The 
method has been recently applied to derive the effective Hamiltonian of the 
single qubit in Ref. [22] ■ The turning point in this case is represented by the way 
in which the two qubits are put into connection. Since the two dots composing 
the qubit are asymmetric, it follows that there are different possible config¬ 
urations. For the purpose of our analysis we consider the two configurations 
shown in Fig. [Q 



Fig. 1 Schematic of the two configurations for the couple of interacting double dot exhange 
only spin qubits. Dotted lines indicate the main interactions. 


Generally speaking the Hamiltonian model is always expressed as the sum 
of three contributions: 

H = H a + H b + H ab . (1) 

The first two terms, common to all the cases under study, are the free Hamil¬ 
tonians of the two qubits hereafter called a and b, written in terms of the 
creation and annihilation fermionic operators c k and Ck respectively 


H q = H eq + H tq + H Uq + H Jq 


( 2 ) 
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where q = a,b and with 


3g 

H eq = e fc c io- c fccr 

fc=lg ,cr 

= *1,3, ^(4,3 C 3,3 + *1-C-) + *2,3, Xi(4,3 C 3^ + /‘■CO 


H u.= H t^fcrzfctrifcj, + !7i,2 5 (ni 5 t + ^i 5 i)(«2,t + ”2,4)+ 


fe=i„ 


+ U lq3q (n lqt + n lqi )(n 3qt + n 3ql ) + U 2q3q (n 2qt + n 2ql )(n 3qt + 713 , 4 .) 


H. 7 . = H 


(Iq3g) ! jj(2q 3 q) ( ^(l g 2g) 


in 


where q = a, b and 

+ Jp <ZJ ' , ^ C jqt C j ? 4. Ci qt Ci q4. + ^ ^ f, CT C jg CT + ^- C -) 


( 3 ) 


( 4 ) 


for every pair of spins considered. H eq and ift describe respectively the single 
electron energy level of each dot and the tunneling energy. The last two terms 
Hu and Hj q constitute the intra-dot and inter-dot Coulomb interactions. 
More in details, the parameters appearing in Eq.© are: the spin exchange 
ji* 9 , the pair-hopping Jp q '’ q ^ and the occupation-modulated hopping terms 

ji^-qjq) 

J t 

The last term appearing in Eq. © is the interaction Hamiltonian H ab , 
whose explicit form depends case by case by the interaction mechanism, but 
always expressed as: 

H ab = H t + H v + Hj. (5) 

In the next subsections we are going to present the final form of the effective 
Hamiltonians in the two cases depicted in Fig. [L] 


2.1 Configuration A 


The configuration A in Fig. [T] corresponds to the situation in which the right 
dot of the first qubit is put into direct connection with the left dot of the 
second qubit in such a way that the principal mechanism of interaction is 
between one orbital of one dot with two orbitals of the other. The interaction 
Hamiltonian m is given by the sum of the following terms 


Hi =t 3a l b 51 ( 4 a<T C Uo + h.C .) + t 3a2b J2( C 3aCT C 2 b a + h.C.) 


Hu =U 3a i b {n 3at + n 3al )(n lb ^ + n lbl ) + U 3a2b (n 3at + n 3al )(n 2bt + n 2bi ) 


Hj —Hj 


(3alb) 


H 


(3a2 b ) 
J ’ 


( 6 ) 
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where Hj is defined as in Eq. Hi- 

Following the same procedure reported in Ref. \T2 \, we calculate the pro¬ 
jected Hamiltonian. The effective Hamiltonian, concerning the low energy ex¬ 
citation, appears as the sum af all the exchange interactions between pairs of 
spins 


H eff = Y {Ji q \Si q ■ S 3 , + J 2 ,3 9 S 2 , • S 3 , + Ji q2q S lq ■ S 2q ) + 

q=a,b 

+ J 3alb S 3a .S lb + J 3a2b S 3a -S 2b . (7) 

The effective coupling constants, considering that intra-dot Coulomb energies 
are larger than all the other contributions are explicitly given in Appendix A. 


2.2 Configuration B 

The second configuration B under study, (see Fig. H, corresponds to the sit¬ 
uation in which the direct interaction is between the dots that present two 
orbitals. This structure create a greater number of interconnections between 
qubits than in configuration A. We assume that energy detuning between the 
double occupied quantum dots is small enough with respect to energy differ¬ 
ence between levels 2 q and l q in each quantum dot. As a result, terms of the 
Hamiltonian containing tunnelling rates and exchange interactions between 
energy levels in different quantum dots with different indexes (i.e. l a 2&, 2 a lb) 
are negligible. The interaction Hamiltonian terms are finally given by 


Ht —ti a i b 


+ h.C.) + t 2a2b Y( C 2 a <T C 2„° + hx -) 


Hu =Ui aU {ni a T + n lai ){n lbt + n lbi ) + U la 2 b (n la1 - + n lai ){n 2bt + n 2bi )+ 
+ U 2 a i b (n 2a i- + n 2 al )(n lb t + n lbi ) + U 2 a 2 b (n 2a t + n 2 ai ){n 2b t + n 2bi ) 


Hj=Hj 


Half, 


H 


(2 a 2 b ) 


( 8 ) 


where Hj is defined as in Eq. 0. 

Analogously to the previous case, the effective Hamiltonian appears as the 
sum of exchange interactions: 


H eff 


Y^ ('^ 1 g 3 9^’ 1 g ‘ + J2 q 3 q ^2 q ' S3, + Jl q2q Sl q ■ S 2() ) + 

q=a,b 


+ ^lalfcSla ‘ + Jl a2fc S la • S 2fc + J 2a l fc S 2a • S l6 + J 2cL 2 b ^2a ' S 2() . 


(9) 


The effective coupling constants, under the assumption of larger intra-dot 
energies with respect to inter-dot ones are explicitly given in Appendix A. 
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3 Gate sequences for exact CNOT gates 

This Section presents the gate sequences which realize CNOT operations in 
both configurations A and B. The effective Hamiltonians, in which only ex¬ 
change interactions between pairs of spin with effective coupling constants 
appear, allow to represent the unitary evolution between the i-th and j-th 
spins with the operator 

U i:j (t) = e~s u « Si ' s h (10) 

After having introduced the mathematical background in Appendix B, we 
examine the simplified toy model in which all the interactions between spins 
are assumed controllable from the external in Appendix C. In particular we 
follow the procedure delineated in Ref. m composed by two steps: in the first 
one we use the central sequence found in Ref. [20] that give a gate operation 
that is locally equivalent to a CNOT operation; in the second step, single 
qubit operations, to be applied before and after the central sequence in order 
to obtain an exact CNOT gate, are found. 

In this Section we solve the realistic problem in which only inter-dot in¬ 
teractions are tunable from the external, while the intra-dot ones are fixed by 
the geometry of the system. To this end we adopt a search algorithm with 
a variable number of steps, showing how it is possible to obtain directly the 
exact CNOT for the configurations A and B. 


3.1 Physical fixed intra-dot interaction 

Let’s now concentrate on the realistic physical situation in which the unavoid¬ 
able intra-dot interactions are fixed and not tunable from the external. In fact 
J 12 does not depend on the tunnelling rates (see Eas. fill and fl6l) whose values 
can span several order of magnitudes and strongly control J 13 and J 23 . We 
assumed max(Ji 3 ) = max(J 2 3 ) = J max and we set a realistic constant value 
for J 12 = J max /2 to model the control ineffectiveness. Moreover, the exchange 
interactions Ji 3 (t) and ./ 23 (f) are assumed to have instantaneous turn-on and 
turn-off. In these hypothesis we derive for both configurations the sequences 
of unitary operations generated by exchange interactions as in Eq. m, whose 
product gives the CNOT operation in the total angular momentum basis. 

3.1.1 Configuration A 

The interaction sequence for an exact CNOT operation for the configuration 
A is calculated by using the search algorithm in the case of fixed intra-dot 
interaction with a final objective function value of 0.001 (see Eq. l28l) . The 
resulting sequence is reported in Tab[I] [23] . After the CNOT sequence, the 
resulting transformation matrix which couples initial and final quantum states 
in the 9x9 subspace is shown in Fig(2J The 4x4 block in the up-left corner 
corresponds to the graphical representation of the CNOT matrix reported in 
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Ea. (l?5l) . Same correspondence is obtained for the 4x4 block in the 5x5 
subspace (not shown). 


Table 1 Gate sequence implementing an exact CNOT gate for the configuration A with 
fixed Ji2=J rnax /2 in both qubits. Interactions in the same time step can be turned on 
together. The “wait” interaction represents only the fixed interactions aia .2 and 6162 with 
no other interactions on. Times are in unit of h/J rnax . 


Step 

Interaction 

Time 

Step 

Interaction 

Time 

1 

wait 

0.500 

17 

a 3 b 3 

0.299 

2 

b 163 

0.866 

18 

<21 <23 

0.136 

3 

a 3 b\ 

1.409 

19 

a2«3 

0.136 

4 

a 3 b 2 

0.550 

20 

a 3 b 1 

0.543 

5 

aia3 

0.194 

21 

wait 

0.907 

6 

U2U3 

0.194 

22 

b 3 b 3 

0.598 

7 

wait 

0.894 

22 

<12 <23 

0.377 

8 

03^2 

0.474 

23 

aia3 

0.377 

9 

< 21(23 

0.033 

23 

b2b 3 

0.736 

10 

CL2CL3 

0.279 

24 

<2362 

0.248 

11 

<21(23 

0.246 

25 

wait 

0.139 

12 

<23^1 

0.433 

26 

b 3 b 3 

0.215 

13 

<13 &2 

0.325 

27 

wait 

0.254 

14 

wait 

0.210 

28 

a 3 bi 

0.014 

15 

«3 fe 

0.139 

29 

<2362 

0.005 

16 

<22<23 

1.155 








Fig. 2 Graphical representation of modulus and phase of the final transformation matrix 
for the CNOT gate in configuration A with gate sequence reported in Tab. [I] 
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3.1.2 Configuration B 


For the configuration B the interaction sequence is presented in Tabj2]|2I] and 
the final transformation matrix, whose 4x4 block in the up-left corner correctly 
corresponds to a CNOT matrix, is shown in Fig(3] Same correspondence is 
obtained in the 5x5 subspace (not shown). 

Table 2 As in Tab[l]but for the configuration B. 


Step 

Interaction 

Time 

Step 

Interaction 

Time 

i 

wait 

1.200 

15 

wait 

1.604 

2 

a2«3 

0.238 

16 

a2&2 

0.121 

3 

aia,3 

0.238 

17 

wait 

1.312 

3 

bib3 

0.068 

18 

a,2b2 

0.168 

3 

CL2b2 

0.429 

19 

wait 

0.400 

4 

aibi 

0.439 

20 

bib3 

0.066 

4 

^2^3 

0.047 

21 

aibi 

1.994 

5 

felf)3 

0.645 

21 

a,2b2 

0.150 

6 

wait 

1.999 

22 

aibi 

0.268 

7 

02^2 

0.242 

23 

bib 3 

0.846 

8 

wait 

0.064 

24 

wait 

0.036 

9 

blb 3 

0.034 

25 

aifei 

0.193 

10 

aifei 

1.151 

26 

wait 

0.294 

10 

^2^3 

0.006 

27 

^2^3 

0.052 

11 

02^2 

0.021 

27 

ai&i 

0.885 

12 

^2^3 

0.271 

28 

a2^2 

0.154 

12 

ai&i 

0.448 

29 

wait 

0.092 

13 

a2^2 

0.111 

30 

b2b3 

0.459 

14 

62^3 

0.099 

31 

wait 

0.125 



Fig. 3 Graphical representation of modulus and phase of the final transformation matrix 
for the CNOT gate in configuration B with gate sequence reported in Tab. [2] 
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Fig. 4 The single qubit device is modeled as a silicon nanowire (in green) embedded in an 
insulator slab (in yellow). Accumulation gates forming the quantum dots are highlighted in 
red whereas the contacts controlling the inter-quantum dot electrostatic barriers are shown 
in blue. The back gate is shown in gray. 


4 Gate design and performances 

In this Section the design of the single qubit holder is presented and evaluation 
of the CNOT gate times are provided by using a Spin Density Function Theory 
(SDFT)-based simulator. The simulator solves the Kohn-Sham equations in 
the Effective Mass Approximation (EMA) with anisotropic effective masses for 
each couple of valleys of silicon along A crystallographic directions and for both 
spin down and spin up populations [25] . When the eigenstates are obtained, 
the spin density concentrations are calculated and the effective potentials, 
namely the Hartree and the LDA exchange-correlation potentials, are derived 
under the Local Density Approximation (LDA). The total potential is then 
calculated self-consistently by solving the Poisson equation with the applied 
potentials from the external. The simulation ends when the error between the 
potential of the current iteration and that of the previous one is under a given 
tolerance. 

The simulated device, presented in FigH] features realistic geometric sizes 
closer to the device (single and double quantum dots) proposed in Refs. [1251 
[27]. The device considered is a silicon nanowire featuring a rectangular section 
with a fixed thickness T$i =15 nrn and with a width W. An AI 2 O 3 layer with 
thickness Tai 2 o 3 =40 nrn is deposed on the nanowire and 30 nrn-wide Al gates 
placed orthogonally to the nanowire direction and separated by di n t. er Gate are 
used to electrostatically confine electrons and control the inter-dot tunneling 
rates in the underneath silicon. 

The design process is the following: first of all the energy level of singlet 
and triplet states of the double occupied quantum dot are calculated and the 
singlet-triplet energy difference AE$t is derived for two different values of W, 
30 and 60 nm, as shown in Fig. [5] AE$t has to be greater than the thermal 
energy in order to avoid unwanted transitions from singlet to triplet states. 

Then to exploit the two lowest quantum states from a single A valley, the 
valley splitting AE V needs to be higher than the AE$t- AE v is defined as 
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Fig. 5 Singlet-t 



Fig. 6 Valley splitting AE V as a function of the back gate electric field Fb ac k for VF=30, 
60 nm. 


the difference between the two lowest energy values of fundamental couple of 
A valleys. AE V is enhanced by increasing the electric field at the Si/A^Oa, 
interface that can be indeed obtained by polarizing negatively the back gate. 
In Fig. [G] the valley splitting AE V is reported as a function of the electric field 
at the back interface Ft ac k- 

Subsequently, the evaluation of the tunneling rate TR as a function of the 
inter-gate distance dinterGate is necessary to evaluate the obtainable exchange 
interaction values and the corresponding gate sequence time t. In Fig. 0 the 
tunneling rate TR is reported as function of dinterGate, showing that a linear 
reduction of dinterGate can increase TR roughly exponentially. 

Summing up the step times of the sequence (see Tables |T] and [2]) gives 
the total gate time in units of h/J max . Given the tunnelling rate TR and 
the singlet-triplet energy splitting AE$t, the maximum exchange interaction 
jmax j g estimated by using the formula J max =TR 2 /AE st ■ As a result, an 
estimate of the total gate time can be calculated. Simulations show that when 
dinterGate= 45 nm, W =60 nnr and Fb ac fc=0.45 MV/cm, the CNOT sequence 
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Fig. 7 Tunneling rate TR as a function of the inter-gate distance di n t er Q a t e with W= 60 
nm. 

durations are t.Q NOT =7A ns and t^, NOT = 10.9 ns for the configuration A and 
B, respectively. From Ref. m a preliminary experimental dephasing time T 2 * 
longer than 20 ns can be extracted, suggesting that at least a CNOT operation 
should be performed successfully. 


5 Conclusions 

A pair of interacting double dot exchange only qubits is exploited to realize 
CNOT operations. We express the evolution operator in a simple form with an 
exchange-only effective Hamiltonian model of the composite system to obtain 
sequences with appropriate interactions and times. The search for sequences in 
two different configurations was performed numerically developing and using 
a mixed simplex and genetic algorithm. We compare the non physical case 
in which all the interactions are controllable from the external and the re¬ 
alistic condition in which intra-dot interactions are fixed by the geometry of 
the system. Imposing realistic conditions to the interactions produces CNOT 
sequences with a number of steps up to 31. 

We compare qubit designs featuring different sizes by exploiting a SDFT- 
based simulator. Gate performances are calculated, providing CNOT gate 
times of t^ NOT =7A ns and tQ NOT = 10.9 ns for the configuration A and B, 
respectively. 


A Details of the calculation of the effective Hamiltonians for the 
couple of interacting qubits 

In this Appendix we are going to report all the detailed expressions for the exchange coupling 
constants between pair of spins in both the configurations examined. The first (last) three 
indices inside parenthesis, 0 < i / j / fc < 2, assuming only integer values, denote the 
number of electrons in each level for qubit a (b). 
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A.l Configuration A 

The coupling constants for the configuration A are given by 


J l q 3 g 

■h q 3 q 

J l q 2 q 

J 3al b 

J 3 a 3 b 


-r^m q 3 q - 4 lq 3 q) ) a - 2 ji lq3 * ) 

A&l q 

^ A(+i r( 2 q3(j)s2 o r(2q3g) 

aeT ^ 9 9 - t e 


1 


1 


AE 3 AE 4 


4j( 1 q 2 g) 2 _ 2 J^ lg2q ) 


1 


4(t 3o i 6 - J t ( 3 alb) ) 2 -2j; 


j(3a Ifo) 


^{t 3 a 2 b -J^) 2 - 2 J^ 2 »\ 


where the energy differences for qubits a and b are 


^^a( b ) = -^( 012 ,in) (-^( 111 , 012 )) ~ ®(lil,lll) 
AE 2 a (b) = -^( 102 , 111 ) (-£'( 111 , 102 )) - 

^ E3 a(V) = ■^'(201,111) (■^'(111,201)) — ^(lii.in) 
^® 4 a(b) = -^(021,111) (■^'(111,021)) — 

AE 5 = E( 112 i oh) — E( ln ln ) 

AEg = E( 112101 ) — E( ln nl ) 


( 11 ) 


( 12 ) 


with 


E (ijk, 111) —t£l a ~h JE2 a + k£3 a ijUl a 2 a + ifc?7l a 3 a + kjU2 a Z a + $i2Ul a + Sj2U2 a + ^k2^3 a ^~ 
+ £ lfc + £ 2 b + £3 6 T^liA + Hl63b + U 2 b 3 b + 

+ kU 3a i b + kUs a 2 b (13) 


E (lll,ijk) — £la + £2 a + £3a + Hl a 2 a + Hla3a + H2 a 3 a + 

+ *£lb + i £ 2 b +^£3 6 +tjtll626 +ifcHl636 + kjU2 b 3 b + Si2U\ b + 5j 2 U2 b + 8 k2 U 3b + 
+ iU 3 a l b + iU 3 a 2 b (14) 


1^(111,111) — £ la + £ 2 a + £3a + Hla2a + E l a 3 a + U2 a 3 a + £ l b + £ 2 b + e 3 b + Ul b 2 b + U\ b 3 b + U 2b 3 b + 
+ U 3a l b +U 3a 2 b (15) 
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A.2 Configuration B 

The coupling constants for the configuration B are given by 

a, 3 , ^ ■ J ‘ (M,))2 “ 2 J e lq3q) 

J2 q3q ^ -7^4(t 29 3, - jT q3q) ? - 2 Ji 2q3q) 

—L_ + 4 J 1 (1,2,)2 - 2J ( e lq2q) 

AE 3q AE 4g ) t 

Ji a i b ^ -2 j( lAlB) 

A a 2 6 = 0 
J 2al b = 0 

J2 a 2 b — —2J^ a2b \ 

where the energy differences for qubits a and b are 

^ El a(b) = £( 012 , 111 ) (£( 111 , 012 )) “ £( 111 , 111 ) 

= £( 102 , 111 )(£( 111 , 102 )) — -£( 111 , 111 ) 

= £( 201 , 111 ) (£( 111 , 201 ) ) “ £( 111 , 111 ) 
^£ 4 a (!>) = £( 021 , 111 ) (£( 111 , 021 ) ) — £( 111 , 111 ) 


(16) 


(17) 


with 

£(ij'fc,lll) =* £ la + jS2 a + ke 3a +ijUl a 2 a + ikUl a z a + kjU2 a 3 a + $i 2 Ul a + 8 j2U2 a +<5fc2£3 a - 
+ £l b + £2 b + £3 b + U\ b 2 b + U lb3b + U 2b 3 b + 

+ *£lalb + *£la2(, + jU2 a l b + jU2 a 2 b (18) 

£(111 ,ijk) =£ 1 a + £ 2 a + £3 a + U la2a + Ul a3a + U2 a 3a + 

+ i£i b +j£ 2 b +ke 3h + ijU\ b 2 b +ikUi b3b + kjU2 h 3 b + Si2Ui b + Sj2U2 b + 5 k2 U 3b 
+ *£lalb +i£la2„ +iU2 a l b +jU2 a 2 b (19) 


£(111,111) —£1 a + £2a + £3a + Ul a 2 a + £l a 3 a + £2 a 3 a + £l b + £2t + £3 b + Ul b 2 b + Ul b 3 b + U2 b 3 b + 

+ U\ a \ b +U\ a 2 b +U2 a l h +U2 a 2 b (20) 


B Mathematical background 

In this Appendix the fundamental mathematical tools used in the following to derive the 
sequences that realize the CNOT gates are presented. 

Let’s introduce the logical basis {10), 11)} used hereafter for each qubit. It is composed 
by singlet and triplet states of a pair of spins, for example the pair in the left dot, in 
combination with the angular momentum of the third spin, localized in the right dot. This 
means that the logical states are finally expressed in this way 

|0> = |So>| \r), |1> = \j\\T 0 )\ i) - y||T_>| t) 


( 21 ) 
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where |<So), |To) and |T±) are respectively the singlet and triplet states, whose explicit form, 
in terms of the eigenstates of cr z , is here reported for completeness 

|So) = 1 U) J W , |T 0 ) = 1 W J, 1 W , |T_) = |U>, |T+) = |tt>- (22) 

The basis state of the Hilbert space containing three electron spins, representing one 
qubit, written in the computational basis via Clebsch-Gordan coefficients is given by: 


|l> = |So>lt> 

|2) = |5 0 >U) 

|3) = -h(V2|T+)U)-|3b)|t>) 

|4) = -^(|T 0 )U)-V / 2|T_)|t>) 

|5> = |T+)| t) 

|6> = T(| T+)U ) + v / 2|T 0 >|t ) ) 

|7> = T (v^|T 0 >U> + |T_>|t>) 

|8> = |T_>| 4.) (23) 


with |So), |To) and |T±) defined in Eqs. (122() . On the other hand the composite system of 
two qubits, that is six electron spins, is hereafter described by a nine-dimensional basis in 
the subspace with total angular momentum operator equal to S = 1, S z = — 1 obtained 
composing the one qubit states in Ea. (l23t) with appropriate Clebsch-Gordan coefficients: 

l^ 9) > = |2> |2> 

l4 9) > = I2)|4> 
l4 9) ) = |4)|2) 
l4 9) > = |4>|4> 

i4 9) > = ^ii)is)-^|2>|7) 

!4 9) > = ^|3>|8>-i|4)|7) 

i4 9) > = -^|8>ii>+^|7>|2> 

l4 9) ) = -^|8)|3> + i|7>|4> 

i4 9)> = \\[\ (|6>|S> + 18>16)) - 17)|7). (24) 

Basis vectors ] ) are valid encoded states and correspond exactly to the logical 

basis {100), 101), 110), 111)} in which the CNOT gate has the usual form 


CNOT = 


/l 0 0 0\ 
0 10 0 
0 0 0 1 
\0 0 1 0 / 


(25) 


basis vectors - |5g 9 ^) are leaked states. 
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The objective function used in the fixed-step genetic algorithm to derive single qubit 
operation is: 

fcNOT = \/l ~ \VIZ) + U (2,2) + U (3A) + ( 2 6) 

where are the matrix elements of the CNOT in the encoded states in the 9x9 subspace. 

The objective function is exactly equal to zero when all the U^) entering into Eq. d 26 D have 
modulus 1 and a common phase in each subspace. 

To study the realistic situation in which intra-dot interactions are fixed by the geometry 
of the system, similarly to Ref. m a search algorithm with a variable number of time steps 
is developed and used. Following the procedure in Ref. m we adopt an objective function 
for the genetic algorithm that, due to the structures of exchange matrices, is confined into 
two subspaces for the total angular momentum operator of the composite system. The basis 
of the first 5x5 block with total angular momentum S = 0, S z = 0 is given by: 

| 6 ®> = -^(| 1 >| 2 >-| 2 >| 1 » 

|4 5) > = T (11)14) -|2>|3» 

|4 5) > = -L(|3)|2>-|4>|1» 

|&i 5) > = -L(|3>|4>-|4>|3» 

l4 5) > = \ (|5>|8> - 18)15) + 17)16) - |6)|7» , (27) 

where the states on the right of the equations are given in Eq. (f23ll : the basis of the second 
block of dimension 9x9 with S = 1, S z = — 1 has been previously defined in Eq. ED- 
The three 9x9 blocks in correspondence to S' = 1 are identical, so they need not being 
constrained separately. Moreover the CNOT matrix on the subspaces in correspondence to 
S = 2 and S = 3 is completely unconstrained, which is automatically satisfied with the 
exchange gates. The objective function is finally defined by 


fcNOT = \Jz- \\+ u% 2) + if A) + u$ 3) \ - i|t/W J + u$ 2) + C/ ( ( 3 9) 4) + C7 ( W } |, 

(28) 

having the same meaning of the objective function in Eq. J26D . In the case of the CNOT 
operation the subspace with S = 0 and S = 1 have the same global phase. 


C Toy model with all controllable interactions 

In this Appendix the CNOT sequence for the toy model with all controllable interactions is 
analyzed. In the first part the uncorrectness of sequence reported in Ref. is demonstrated, 
while in the second the sequence starting from our genetic algorithm is derived. 

In the following we present the sequence of exchange operations as done in Ref.l20| for 
configuration B, in which a locally equivalent CNOT gate is presented equipped with single 
qubit operations to obtain an exact CNOT. In Fig[8|we represent graphically the modulus 
and phase of the final transformation matrix presented in the supplemental material of 
Ref. 1201 for the configuration B described by the Hamiltonian in Eq. ((9j). It includes the 
sequence to obtain a locally equivalent CNOT and single qubit operations to obtain exactly 
a CNOT gate. 

While the locally equivalent CNOT presented in the supplemental material of Ref. m 
returns (0,1) Makhlin coefficients m as it should be, when the single qubit operations 
proposed are applied before and after the central sequence, the 4x4 block, contrarily to 
they claim, does not represented an exact CNOT as Fig [8] witnesses. 
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101 

00> 


Fig. 8 Graphical representation of modulus and phase of the complete transformation 
erroneously called exact CNOT gate in the supplemental material of Ref. W . 


Starting from the central gate sequence of Ref. m for the local equivalent CNOT in 
configuration B and by adopting a genetic algorithm with a fixed number of time steps to 
find the single qubit operations with the objective function (12611 . we obtain the exact CNOT 
shown in Fig. [9] 


phase 



Fig. 9 Graphical representation of modulus and phase of the final transformation matrix 
for the exact CNOT gate in the toy model. 


The sequence of single qubit interactions considered is reported in Tab. [3] 
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Table 3 Single qubit operations transforming the locally equivalent CNOT (LECNOT) to 
an exact CNOT. ai, ( 12,(13 (bi,&2 5 &3) denote the three spins in the qubit a (b). Interaction 
times are in units of Sequences on the left and right must be applied before and 

after the LECNOT sequence, respectively. 


Interactions 
before LECNOT 

Time 

Interactions 
after LECNOT 

Time 

<2203 

0.2784 

«2<23 

0.3319 

^2^3 

0.9965 

&2&3 

0.8270 

a\Cl2 

0.4733 

CL\CL2 

0.5270 

t>l&2 

0.2948 

&1&2 

0.4127 

0203 

0.6687 

<22(23 

0.7221 

^2^3 

0.5976 

&2&3 

0.6380 
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